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Abstract 



The method of Morse theory is used to analyze the distributions of unit 
charges interacting through a repulsive force and constrained to move on the 
surface of a sphere - the Thomson problem. We find that, due to topological 
reasons, the system may organize itself in the form of pentagonal structures. 
This gives a qualitative account for the interesting "pentagonal buttons" dis- 
covered in recent numerical work. 
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I. INTRODUCTION 



The Thomson problem consists in determining the minimum energy configurations for a 
given collection of N unit like-sign electric charges constrained to move on the surface of a 
two-sphere of radius R [1]. The Coulomb energy of such a system is 

where fi is the radial coordinate of the i^^ charge on the sphere. The actual values of the 
electric charge and of the dielectric constant of the medium are inessential for the qualitative 
considerations we are going to make throughout the paper. 

This is an old and largely unsolved problem that, in its generalized version - i.e. for 
more general repulsive potentials as well as for topological defects rather than unit electric 
charges - finds applications ranging from superconductivity to biology. 

In recent years, triggered by the seminal paper of reference [2], there has been a lot 
of interesting theoretical, numerical and experimental work on the generalized Thomson 
problem [3]. In [2] the authors deal with disclination defects constrained to move on a 
sphere. They show how an effective defect model, rather than the analysis of the actual 
elementary charges interaction, proves to be extremely reliable to describe (by numerical 
means) the ground state configurations in terms of the ratio of the defect core energies to 
the Young modulus. This way they provide interesting empirical solutions to the Thomson 
problem, presenting patterns the defects form on the surface of the sphere at the various 
energy thresholds. 

We would like to understand the patterns formed the way described in [2], i.e. why these 
charges (defects) arrange themselves on the surface of the sphere following precise symmetry 
prescriptions. We would expect the distribution to have the property that each charge has 
the same environment. If such an assumption is made the charges should be distributed in 
a way that represents a tiling of the surface of the sphere. 

The general problem of tiling a genus g surface, say S^, can be solved in a simple way. 
Suppose we want to tile with regular ]?-sided polygons (p-gons) assembled in such a way 
that each vertex is shared by (is a common vertex for) 3 p-gons, and each edge is shared by 
(is a common edge for) 2 p-gons. If kp is the number of p-gons used^ to tile S^, the resulting 
polyhedron P has Vp = l/SJ2k kpP vertices, Ep = 1/2 X^a: Kp edges, and Fp = J2k K faces, 
giving for the Euler characteristic xi^g) = Vp — Ep + Fp, the following expression 

^ kp{6 -p) = 6x(S,) = 6(2 - 2g) . (1.2) 

k 

Generalization to the case of vertices shared by not always three p-gons can be easily con- 
structed. 

For the sphere = 0, and if only hexagons and pentagons are used, there can be an 
arbitrary number of hexagons, but there must be exactly 12 pentagons. This follows from 
Eq.(1.2) 



There are, for instance, 5-gons (pentagons), k^ 6-gons (hexagons), k-j 7-gons (heptagons), etc. 
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Yl - p) 



h (6 - 5) + kQ (6 - 6) = 12 



(1.3) 



k 



hence = 12 and fcg is arbitrary. If one also uses heptagons for the tihng of the sphere, 
then one has the interesting result 



Thus, one can start off by tiling the sphere with an arbitrary number of hexagons and exactly 
twelve pentagons. Then one can go on by adding an arbitrary number of pairs pentagon- 
heptagon, but not a pentagon or a heptagon separately. Let us notice en passant that other 
ways of tiling the sphere are, of course, possible. For instance, if only equilateral triangles 
(3-gons) are used, then from Eq.(1.2) follows 



that is the tetrahedron. 

We are particularly interested in the 5-gon-6-gon-7-gon tiling because the authors of [2] 
consider the pentagons and heptagons as disclination defects in a sea of hexagons. They 
construct an effective theory for Eq in (I.l) in terms of interactions between the defects, 
the background charge distribution of hexagons merely providing the value of the effective 
Young modulus. They prove that for a sphere of large radius, distorting the curvature 
by the introduction of defects is not excessively expensive in terms of bending energy. A 
5-gon defect makes the local curvature negative, while a 7-gon makes it positive. For a 
large spherical system, therefore, the introduction of a 5-gon-7-gon pair is an energetically 
reasonable way of joining a group of 12 charges. A calculation had to be carried out to see if 
such a charge configuration was energetically favored as compared to 12 charges organized 
on two additional hexagons. 

They find that when the total number of charges on the sphere exceeds a certain critical 
value of (9(500), the system prefers to have a collection of 5-gons and 7-gons arranged in the 
form of a linear chain of alternating 5-gon-7-gon sequence ("scar"). For even larger values 
of A^, the preferred form is what they name "pentagonal buttons" . These are configurations 
of two nested circles, with five 5-gons placed on the outer circle, five 7-gons on the inner 
circle, and finally a 5-gon in the common center. As N further increases, the defect system 
formes more intricate patterns. A Cs-symmetric configuration of defects is noted. 

In this paper we would like to give qualitative topological arguments suggesting why these 
different defect configurations - such as pentagonal buttons or C3 symmetric configurations, 
discovered by numerical minimization of Ec - can appear. We shall do so by recalling in 
Section II how the topological methods of Morse [4] apply to second order phase transitions 
in crystals (as pioneered by Michel, see e.g. [5]) paving the way to the application to the 
Thomson problem we shall deal with in Section III. Section IV is devoted our conclusions. 

II. SYMMETRY BREAKING, MORSE THEORY AND SELECTION RULES 

A crystal has a density function p{x) invariant under one of the finite crystallographic 
groups, say G 



h-kT^l2 . 



(1.4) 



A;3 3 = 12 or A;3 = 4 , 



(1.5) 
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p{gx) = p{x) with g eG , 



(11.1) 



where p{x)dV is the probabihty to find an atom of the crystal in the volume dV. Not all 
finite groups acting on 3 spatial dimensions, x — (xi, ^2, X3), are crystallographic groups, 
i.e. actual symmetries of crystals in nature. For instance: Oh, the octahedral group, the 
largest symmetry group of the cube, is a crystallographic group, while /, the icosahedral 
group, largely used in [2] for the effective theory of Ec, is not. 

Crystals undergo second order phase transitions. The theory of such transitions has been 
estabhshed by Landau [6] . The density function p changes smoothly from one phase to the 
other, while the symmetry group G suddenly changes to a subgroup H <Z G. The density 
function can then be decomposed as 

p{x) + 5p{x) {11.2) 

where p is G-symmetric, while 5p is iJ-symmetric, and, at temperatures below the critical 
value T < Tc, 5p — while for T > Tc is (^p 7^ 0. As the function changes continuously 
during the transition, 5p is small near Tc- 

These transitions are (proto-)typical examples of spontaneous symmetry breaking: The 
thermodynamic potential ^{p,T,P) is always G-invariant^. For T < Tc the vacuum is G- 
invariant, realizing an explicitly symmetric phase; for T > the vacuum is iJ-invariant, 
realizing a spontaneously broken symmetric phase. By vacuum here we mean the p config- 
uration that minimizes the functional $(p) for the given values of T and P. 

Thus, given a G-symmetric crystal (p and $) one can select the sub-group H into which 
the crystal will make the transition, among the allowed subgroups of G, by finding the 
minimum of $ with respect to p. Using physical reasonings Landau assumed that, near Tc, 
$ is a polynomial in p of three terms of order zero, two and four (the Landau polynomial). 

As p is a G-invariant function, it can be expanded in the basis of the functions 
{0j(a;)}i=i,...,ord(G)) where ord(G) is the number of elements (order) of the finite group G, 
and 

ord(G) 

US) = E D,,[g]<PAS) ^ = 1> •••,ord(G) , (II.3) 

for any g E G. The matrices {Aj} reahze a reducible representation of G. Being p real, the 
representation is necessarily orthogonal. 
In particular 

ord(G) 

Sp{x) = J2 ^r^S) 
i=l 

and 

ord(G) 

5p'{x) = 5p{gx) = E V'M^) (n.5) 
1=1 



^This is so because it cannot depend on coordinates, hence, in particular, is invariant under G. 
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with ry'j = J2j Dij[g]'i]j. The rjiS are the order parameters of the phase transition {rji = 
at T < T,.. rji 7^ at T > Tc), hence the thermodynamic potential will be minimized with 
respect to these quantities $(r/). The reducible set {rji} can be decomposed into irreducible 
sub-sets where 

n 

vt-T^daMVi a^l,...,n, (II.6) 

6=1 

the dab[t/]s being the n x n irreducible matrix-blocks in the ord(G') x ord(G) matrices Dij[g] 
and, supposing there are m such irreducible representations of dimension ni,n2, ■■■,nm, we 
have 

nl + nl + ... + nl^ord{G) . (II.7) 

Note that {rj^} (or the corresponding {0"}) is not a complete set. Nonetheless, it is only 
one such sets that is needed to expand 6p near Tc [6]. Thus, the problem of finding which 
way the G-symmetry of the crystal is spontaneously broken down to the smaller symmetry 
of one of its subgroups - say H - boils down to the minimization of the real function(al) 

^■.ffeS"^^ff}eR (11.8) 

where n is one of the m values in (II. 7), and we added to i?" the point at infinity that, for 
stabiUty, has to be included as a maximum, so that i?" + {00} = S"- . In what follows we 
shall take n — 3. 

One way of solving this problem is, of course, to explicitly know the coefficients of the 
Landau polynomial that depend upon the details of the model. Many things, though, can 
be said about the critical points (minima, maxima, saddle points) of a real function like 
$ in (II. 8), without knowing its explicit expression, but rather exploiting the constraints 
associated with the topology of yS", as shown in the mathematical work of Morse [4] (for a 
physically intuitive introduction see, e.g., [7]; for extensive applications to phase transitions 
in crystals see, e.g., [5]). 

In a nutshell, Morse proved that, when a function is like $, i.e. smooth, real and defined 
over a compact differentiable manifold (as said, we shall focus on the case of S^), then the 
following constraints hold 

Co > bo (II.9) 
ci-co>6i-6o (11.10) 
C2 - Ci + Co > 62 - &i + bo (11.11) 
C3 - C2 ci - Co = 63 - ^2 + ^1 - bo (11-12) 

where c^ is the number of critical points of $ with index^ i [i = are the minima, i = 3 the 
maxima, < £ < 3 the saddle points), and the 6^s are the Betti numbers of S^. Recall that 



^In the neighborhood oi xq, a result of Morse theory states that the function / has a local 
description of the form 
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bo{S') = 1 = bs{S') , (11.13) 

while the other Betti numbers of are all zero. This is a powerful result: without knowing 
the actual form of the function, the topology of the manifold tells us that for each £ there 
are at least bi critical points with index £ 

ci>bt. (11.14) 

For the case in point we know that $ has at least one minimum (6o = 1) and at least one 
maximum (&3 = 1). Furthermore, denoting by ff^^ e the critical point of $ with index £, 
and by iJ^^^ the subgroup of G that leaves it invariant, 

^f = E46[%f /lei/W, (11.15) 

6=1 

the number of elements (order) of H^^^ is simply related to q 

ord(i/(^)) = ord(G)/Q . (11.16) 

It is a result of the theory of finite groups that ord(i/W) as written in (11.16) is indeed an 
integer, the proof being based on the coset decomposition of G with respect to iJ*-^-* [8], [7]. 

The subgroup we are looking for is the symmetry of the vacuum configuration ff^^ 
and new symmetry H of the crystal beyond Tc- The constraints (II.9)-(II.12) and (11.16), 
the requirement to have at least two maxima (at ff^^ = oo and at ff^^ = 0), i.e. C3 > 2, and 
the fact that the number of real solutions of d^/df] — (with f) & and $ a fourth order 
polynomial) is bounded to be 

Co + ci + C2 + ca < 3^ + 1 = 28 , (11.17) 

prove to be enough to identify the permitted if^°^s! Which one is the actual H of the crystal 
after the spontaneous breaking of symmetry G ^ if, of course, depends upon the particular 
physical situation under investigation. 

Examples are crystalline structures, such as the alloy Cu-Zn, that exhibit the octahedral 
symmetry, whose group Oh has order 48. The subgroups of Oh allowed by the requirement to 
be crystallographic groups, have order 8, 6, 4 and 2. Carrying on the analysis outlined above 
furnishes: OTd{H^^^) = 6, which corresponds to C^v, and OTd{H^^^) = 8, which corresponds 
to C/i,; (recall that ord(Cn) = n and ord(C„t,) = 2n). This means that C^^ and are 
selected as the only possible symmetries of a crystal, originally O/j-symmetric, undergoing a 
second order phase transition. More on this can be found in [5] and in [7]. 



fix) = f{xo) -yl-yl y| + yf+i + • ■ ■ + ydimM > 

where yi = {x — xo)i is the i^^ coordinate of the vector x — xq, and dimM is the dimension of the 
manifold M. The index £ of the critical point xq refers to the number of negative terms present in 
the local description of /(.x). Then £ = means that no negative term is present in the local 
representation for /, so that the critical point is a local minimum. Similarly if ^ = dimM, then 
it means that / decreases in value in all the dimM directions of the manifold, hence this critical 
point is a maximum,. 
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III. THE THOMSON PROBLEM IN THE ICOSAHEDRAL APPROXIMATION 



We shall now extensively rely on the model built up in [2], which deals with a system 
of topological defects on the surface of a spherical elastic material. Our considerations will 
be qualitative in nature, hence we shall skip many of the details of that analysis and will 
focus only on the topological and group-theoretical aspects of the problem. This way we 
shall provide an explanation of (some) of the empirical results of that work with the hope 
that our analysis could provide useful insight for that problem and, more generally, for the 
Thomson problem. 

The authors of [2] deal with a sphere, originally tiled with an arbitrary number of 
hexagons and exactly twelve pentagons. The topology of the sphere requires that the defects, 
emerging as disclinations in the underlying structure of hexagons, must be: first the initial 
12 pentagons, then the pairs pentagons-heptagons. On this we commented in the Introduc- 
tion. For interactions among the defects that are strong compared to the interactions among 
the particles forming the material, as collectively described by the Young modulus, (what in 
[2] is identified as the "large core" energy regime), the twelve pentagons are the only defects 
in the system and arc at the vertices of an icosahedron, because this configuration minimizes 
the energy^ Ec- In the "small core" energy regimes new defects are energetically allowed 
as 12-plets of pentagon-heptagon pairs, and again they prefer to arrange themselves at the 
vertices of icosahedra. The patterns these additional defects describe on the surface of the 
sphere is our concern. 

To try to understand that, we shall apply the analysis of the previous Section by con- 
sidering the energy as a functional of the density function. Thus we shall regard Ec as a 
map 

Ec-.fjeS'' ^ R, (III.l) 

where, as explained, n is related to the symmetry group, the icosahedral group / in this case, 
the point at infinity is added for stability as a maximum, and the notation is the same as 
before. In the large core energy regime Ec and its minimum^ (vacuum) are /-symmetric. In 
the small core energy regimes we shall make the approximation that Ec is still /-symmetric^, 
but the vacua will only be invariant under a subgroup of /. This way we are dealing with 
the spontaneous breaking of the icosahedral symmetry into smaller symmetries we want to 
determine. 



^We shall still call the energy Ec- This is a clear abuse of notation because the energy here is 
quite more complicated than the expression given in (I.l). Nonetheless, for our considerations the 
actual form of the function is not important. 

^Recall that (in the language of the Morse inequalities) cq = ord(/)/ord(/) = 1. 

^That the /-symmetry of Ec in the small core energy regimes is a good approximation is a result 
we borrow from [2], where it is shown that indeed the energy is not severely modified by the 
introduction of these new defects. 
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The icosahedron is a polyhedron with twelve vertices, twenty faces and thirty edges. 
Its symmetry group I has 60 elements arranged into 5 conjugacy classes (the number of 
elements in the class is in parentheses) [8]: Identity (1), C2 (15), C3 (20), C5 (12), C| (12). 
The subgroups of / are C2, C3, C5 (of order 2, 3, 5, respectively) and those obtained by the 
coset decomposition of / with respect to them. The list of the orders of the subgroups is 

{2,3,5,12,20,30}. (III.2) 

Of course, by definition of polyhedron, the vertices of the icosahedron he on the surface of 
a sphere. Nonetheless, we want to stress this here because, together with the topological 
constraints on the tiling and the choice of the irreducible representation of / relevant for the 
problem in point (i.e. the three-dimensional one), it makes clear the role of the two-sphere 
in configuration space. These comments are also in order because the topology of a different 
sphere, the S'^ — in (III.l), will now play an important role. 

Near the phase transition we suppose that Ec is a polynomial of order four, just as for 
the thermodynamic potential $ of crystals. Hence the number of real solutions of E'c = 
is bounded from above 

Co + ci + C2 + C3 < 3=^ + 1 = 28 , (III.3) 

where the point at infinity is also added. 

We can now use the relevant Morse inequahties for the function in (III.l) with n = 3 

Co > 1 (III.4) 

Co < 1 + ci (III.5) 

Co + C2 > 1 + Ci (111.6) 

Co + C2 = Ci + C3 . (Ill- 7) 

These, the request that C3 > 2 and the constraint (III. 3) are what we need to identify the 
order of the allowed subgroups H^. We just have to insert the given order from the list 
(III.2) into 

Co = 60/ord(i/°) . (III.8) 

The resulting list of allowed subgroups consistent with the Morse inequalities and with the 
constraint (III. 3) is 

{5,12,20,30}. (III.9) 

Let us analyze the case ord(iJ°) = 12. In this case the /-symmetry has been broken in 
such a way that there are 5 distinct minima (co = 5) connected by a C5 transformation 
(recall that 12 is the order of the coset of / with respect to C5). Each minimum corresponds 
to a new icosahedron of defects, thus, according to this argument, there would be 5 new 
icosahedra. But, as shown, the new defects cannot appear unless they are in pairs of 5- 
gons-7-gons, thus, being 5 odd, we must have 10 new icosahedra: 5 icosahedra of pentagons 
next to 5 icosahedra of heptagons. These are the pentagonal buttons found in [2] . It would 
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be interesting to see how the residual symmetry of order twelve is related to the way these 
buttons arrange themselves on the sphere which, according to [2], is at certain vertices of 
the rhombic tricontahedron. 

Similarly the case ord{H^) — 20 refers to 3 distinct minima for the new symmetry, 
connected by a C3 transformation. By means of the same constraint to have 5-gons-7-gons 
pairs, we would expect to see 6 new icosahedra at the time, while the numerical results show 
7 such clusters. This discrepancy could be due to the approximations we are introducing, 
either for the invariance of the energy or of the order of the polynomial near the phase 
transition. 

We find also the cases: ord(if°) = 30, which refers to 2 distinct minima connected by a 
C2 transformation, and the case ord(i?°) = 5 which means H — C5, in the notation of the 
previous Section. 

IV. CONCLUSIONS 

We applied the method of Morse theory and second order phase transition in crystals 
to the Thomson problem in the icosahedral approximation of reference [2]. We are able to 
describe the pentagonal buttons and the C3 symmetry found there, although this last one 
only to an approximated extent. We also notice residual symmetries of orders not seen in 
that numerical work. This does not necessarily means that these new symmetries are to 
be found exactly in form given, because the analysis we made might require further care in 
order to be apphed to other cases. 

Our results are qualitative and preliminary. For instance, within this approach we are not 
able to address the important issue of the various energy thresholds involved with the phase 
transitions. Thus we cannot say when the transitions happen. Nonetheless, we can address 
the question why these transitions occur and in such a general fashion that it is suitable, in 
principle, for a wide range of applications. Hence, we believe, the results presented here are 
interesting especially in the view of stimulating the discussion for a deeper understanding 
of the Thomson problem in general. 

We also notice that, having clarified the topological origin of the pentagonal buttons and 
Ca-symmetric configurations, it is reasonable to argue that these configurations share some 
features of solitons, in the sense that they cannot be undone by a continuous transformation. 
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